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Constraints on Extragalactic Point Source Flux from Diffuse Neutrino Limits 
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We constrain the maximum flux from extragalactic neutrino point sources by using diffuse 
neutrino flux limits. We show that the maximum flux from extragalactic point sources is 
E2(dN^/dE) < 1.4 X 10"'' (L„/2 x lO"*^ erg/s)^/^ GeV cm'^ s"^ from individual point sources with 
average neutrino luminosity per decade, L^. It depends only slightly on factors such as the inho- 
mogeneous matter density distribution in the local universe, the luminosity distribution, and the 
assumed spectral index. The derived constraints are at least one order of magnitude below the cur- 
rent experimental limits from direct searches. Signiflcant constraints are also derived on the number 
density of neutrino sources and on the total neutrino power density. 
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The origin of ultra high energy cosmic rays (UHECR) , 
is still unknown. Active galactic nuclei (AGN), gamma- 
ray bursts (GRBs), or processes beyond the standard 
model have been hypothesized to be the sources of UHE- 
CRs. If nearby AGN are the sources of the highest energy 
cosmic rays and if AGN emit neutrinos in addition 
to photons, protons and other charged particles at com- 
parable fluxes, then individual AGN may be observable 
by current generation of neutrino detectors. However, 
only the nearest sources would be detectable as point 
sources, while the contribution of an ensemble of unre- 
solved extragalactic sources would generate a diffuse flux 
of neutrinos. There are plausible but speculative reasons 
to expect a correlation between sources of cosmic rays 
and sources of neutrinos. Several models predict a dif- 
fuse neutrino flux from AGN, in particular neutrino pro- 
duction has been predicted from the core of radio-quite 
AGN 0,11, and from AGN jets and radio lobes 0,i,[H. 
Direct searches for diffuse fj\ and point flux [sl] by current 
telescopes have set the most stringent upper limits, but 
generally have not reached the sensitivity required, and 
the models suggest that challenges exist even for next 
generation telescopes. One of the primary motivations 
for the construction of neutrino telescopes is to search 
for unexpected sources with no obvious connection to the 
power emitted in the electromagnetic band. 

We show in this paper that the z^-flux from extragalac- 
tic point sources can be constrained by the measured 
diffuse ;y-flux limits, and we also use these results to con- 
strain the neutrino intensity predicted in models from in- 
dividual sources. The derived constraints are one order of 
magnitude below current experimental limits from direct 
searches for energies between TeV-PeV, and below cur- 
rent limits and sensitivities of km^ neutrino telescopes, 
such as IceCube, for energies between PeV-EeV. Since, 
the constraints scale with the power of 2/3 of the mea- 
sured diffuse flux, an expected factor three improvement 



in the diffuse flux sensitivity for 1 year of IceCube Q 
data improves the constraints by another factor two. 

Point sources of neutrinos are observed when several 
neutrinos originate from the same direction, and in the 
context of this study, only the very nearest of an ensemble 
of extragalactic sources are detectable as point sources. 
The number of detectable (or resolvable) point sources, 
Ns, presented in [lo(, is determined for a given diffuse 
v-Omx limit and point source sensitivity. The Ng cal- 
culation is based on three assumptions: (1) the sources 
are extragalactic and uniformly distributed in space; (2) 
the neutrino luminosity follows a power law or broken 
power law distribution; (3) the sources are assumed to 
emit neutrinos with an E~'^ energy spectrum. Later, we 
discuss the robustness of the constraint by investigating 
the validity and caveats of the assumptions. 

The number of resolvable sources for a distribution 
of luminosities Li, per decade in energy is given by: 
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where the parameter ^ which is close to unity, depends on 
cosmology and source evolution as described in [10] . The 
neutrino luminosity of the source, L^, has units of (erg/s), 
and {E,nin, Emax) defines the energy range of the flux 
sensitivity, where Emax = 10'^ i?mm for a typical exper- 
imental condition. For canonical energy spectrum pro- 
portional to E~^, we use the Ultra High Energy (UHE) 
results for all-flavor diffuse flux limits from AMANDA |3| 
to obtain the diffuse i/^-flux: Kdiff = ^^^v^ — (1/3) * 
E2$^_^,, = (1/3) * 8.4 X 10"* GeV cm-^s-^sr"! = 
2.8 X 10^* GeV cm~^s~"'^sr~-'^ valid for the energy interval 
of 1.6 PeV < E < 6.3 EeV. This is the energy interval of 
interest for cosmic ray interaction with energies above the 
ankle. For neutrinos at the Very High Energy (VHE), we 
also use limits from AMANDA [uf, Kdiff < 7.4 x 10"® 
GeV cm-2 s'^ sr-\ valid between 16 TeV to 2.5 PeV. 
So, similar diffuse flux limits exist for the entire inter- 
val from TeV to EeV energies. C point is the experimen- 
tal sensitivity to z/-fiuxes from point sources for an E~^ 
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FIG. 1: Constraints on neutrino point fluxes derived from the UHE diffuse i^-flux limit f?!, and from VHE limit ^T], for two 
representative (L^) = (10""', 2 x 10''^) erg/s. Current AMANDA limit [3], IceCube sensitivity j9|] to neutrino point fluxes 
and IceCube limits [l3l | on fluxes from two individual point sources are also shown (thin solid lines). A sample of model 
predictions for i^^-point flux from extragalactic sources are displayed in thin dotted-dashed lines, which are proportional to an 
spectrum or follow a broken power law. Emission from AGN jet, calculated for a 3C279 flare of 1 day period [3C279] 20|; 
Spectra predicted for IVIkn 501 during the outburst in 1997 [Alkn 501] [2^ and core emission due to pp interactions [3C273] 17|; 
radio-quiet AGN [RQQ] 23]; emission from Cen A as described in \Cen A] [2^, [Cen A] [26|] and [Cen A] [27j; emission from 
M87 [M87] 25], and emission from Coma galaxy cluster [Coma] [2a ]. 



spectrum, where the sensitivity from AMANDA [sl is 



point 



E^(dN^/dE) < 2.5 X 10- 



The diffuse flux Kdi / / parameter and the point flux 
sensitivity Cpomt are linearly correlated by the following 
equation: 



■iivKdiff = 



1 



N, 



X a 



point 



(2) 



where (c/Hq) represents the Hubble distance given by 
c/Ho = 3 X 10^ (km s~^)/77 (km s^^ Mpc~^) - 4 Gpc. 
For the case of iV., < 1 the distance ratio {c/HQ)/rjnax > 
1, which occurs for sources well within the Hubble dis- 



tance. The parameter 



defines the maximum observ- 



able distance for a point source of luminosity L^, which 
is given by: 



47r YlilEjjidx / -^rnin^ Cpoint 



1/2 



(3) 



The constraint also holds for time variable sources, since 
it depends only on the observed luminosity and is inde- 
pendent of the duration of the variability 12] . Similarly, 
it holds for beamed sources, such as GRBs. However for 



luminosities of the order of 10^^ erg/s typical of GRB 
emission, a dedicated search for GRBs leads to more re- 
strictive limits 14]. 

We derive an upper limit on the maximum neutrino 
power density V'^ independently of the number density 
of sources, given by: 



)Kd^ff^3Ax 10 



45 



erg/s 
Gpc3 



(4) 



which is one order of magnitude below the power required 
to generate the energy density of the observed extragalac- 
tic cosmic rays [isj . 

A numerical value for Ns can be estimated by incorpo- 
rating the diffuse j/-flux limit and the sensitivity to point 
sources in Eq. [T] TV, ~ (3.7 • lO-^^cm-i) x {Kd^ff ) x 
{Cpotnt)^^^'^ X (Lagn)^^'^ X 1/^ ~ 10"^ computed assum- 
ing Lagn — 2 X 10"^^ erg/s. We chose to scale the value 
for the neutrino luminosity Lagn = V^/ng = 2 x 10'*'^ 
erg/s for a number density of Ug ^ 10^ Gpc^"^ character- 
istic of AGN sources. The parameter — £,agn — 2.2 ac- 
counts for the effects due to cosmology and source evolu- 
tion that follows AGN [10]. The estimate for ~ 10~^, 
which is compatible with the non-detection of any point 
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FIG. 2: Same flux constraints as Fig. [T] compared to a sample of model predictions for f^-point flux from extragalactic sources 
are displayed in thin dotted-dashed lines, which strongly differ from an spectrum. Emission from 3C273 predicted by 
[3C273] [iq ]. including pp and P7 interactions [3C273] [3]; core emission due to P7 interaction [3C273] [HI; AGN jet continuous 
emission [3C279] [Til; emission from NGC4151 by [AfGC4151] [11] and core emission from NGC4151 due to p-y interaction 
[iVGC4151] [iiSpectra predicted for Mkn 421 [Mkn 421] [3, and blazar flaring Mkn 501 [Mkn 501] [H]; GeV-loud blazars 
[GeV blazar]g^. 



sources. The constraint on i/-flux is determined by re- 
solving at least one source, i.e. by setting = 1 and 
inverting Eq. [T] to solve for Cpoint ■ 
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dE 
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diff 
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1/3 



GeV 



2/3 



(5) 



valid for the same energy range 1.6 PeV < E < 6.3 
EeV of the diffuse flux limit Kdiff- This result de- 
fines a benchmark flux constraint <i>c = E^(dN^/dE) < 
1.4 X 10~^ GeV cm~^ on neutrino fluxes from in- 
dividual extragalactic point sources that produce the 
power required to generate the neutrino flux with L^, = 
2 X 10''^ erg/s. The benchmark flux constraint is 
one order of magnitude lower than present experimental 
limits from direct searches, and strengthen for ensemble 
of sources that generate less power. These results show 



that the likelihood of detecting neutrino signal from AGN 
sources will be a challenge for next generation km-scale 
neutrino telescopes. 

Fig. [1] shows the benchmark constraint on extragalac- 
tic point source fluxes derived from the UHE and VHE 
diffuse flux limits. Models are shown with an energy 
spectrum proportional to E^^ (or approximatelly pro- 
portional over the UHE and VHE energy interval). The 
model predictions can be compared to the derived bench- 
mark constraint, <I>c, by assuming that the specific pre- 
diction characterizes the mean flux ^^°d.ei^ ^^^^^ energy 
distribution from an ensemble of sources. By comput- 
ing the ratio TZ = $c'/$™°'^'=', models are constrained if 
TZ < 1. The results from the constraint <I>c compared to 
a number of models of neutrino point fluxes from extra- 
galactic sources are summarized in Tab. HI 

Models shown in Fig. [2] strongly deviate from an E~^ 
spectrum and in this class of models a direct comparison 
with the benchmark Hux is less straightforward. For 
the models 0, [l^, IH, [131 , the predicted energy spectra 
are integrated over the UHE (VHE) energy interval to 
obtain the total number of neutrinos for the given model. 
The result is compared to the integrated neutrino events 
Nc determined by the benchmark flux <i>c and by the 
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TABLE I; Summary of models for i/^ point fiux from extragalactic sources constrained by the results from this work. The 
benchmark flux <l>e deflnes the flux constraint for an spectrum, which is directly compared to the predicted neutrino flux 
for a given model, The redshift of the source is from Q^, and the parameter ds deflnes the distance of the source in 

Mpc computed according to the relation — z x c/Hq. The neutrino luminosity is computed from (gee text for 

details). Upper bounds on the number density, Us, are given in units of Gpc~^. The ratio TZ — ^cj^^"'^'^^ < 1 determines a 
model constrained by this work. 
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detector neutrino effective area A^f / : 

Nc^tu^e J ^cAeffiE,)dE (6) 

Similarly, the number of neutrino events expected from a 
given model, Nmodei, is computed by substituting the 
predicted energy spectrum for in Eq. [HI The ra- 
tio Nc/N^odei is found to be 0.07, 0.2, 0.03 and 0.17 
for [iVGC4151] [3, [3C273] 0, [GeV hlazar] [H and 
[Cen A\ respectively. 

The maximum number density of extragalactic sources 
rig can be expressed in terms of the neutrino luminosity 
L^, using the relation |^12]: 

ris < 47r ln(--; ) x — (7) 

The number density is inversely proportional to the neu- 
trino luminosity Li, and scales linearly with the measured 
diffuse flux K^iff. Therefore we can set a constraint on 
the number density Us based on the measured diffuse 
flux limits Kdiff, as shown in Fig. [3] The thick solid 
line shows the constraint on cx Kdiff/Li^, so stronger 
diffuse flux limits constrain the neutrino source density 
Us to lower values. The thin parallel lines beneath it 
correspond to improvements in the experimental diffuse 
limit Kdiff by factor of 10 and 100, respectively. The 
experimental sensitivity to point flux Cpomt can also be 
expressed in terms of the number density rig as follows: 

= 3^4^^ f ln§^) x {'^f/^ (8) 

Since the Ug scales as {Lj,)^^/"^ the upper bounds set by 
direct point searches (thick dotted line) have a steeper 



slope compared to the diffuse flux constraints. For neu- 
trino luminosity of bright extragalactic sources with val- 
ues < erg/s, the upper bounds on set by the 
diffuse flux are few orders of magnitude below the bounds 
reached by direct searches. 

We derive limits on the number density for specific 
source predictions if ^™°<^'^'' is assumed to character- 
ize the average flux for an ensemble of similar sources. 
The limits on rig are shown as points in Fig. [3] and 
are summarized in Tab HI The neutrino luminosity per 
energy decade is computed from the source distance 
dg and the flux <j)™o'^'=' using the relation in Eq. [3l 
= $mode' X 4'Kd'j. ln{Emax / E^in) , which assumes 
isotropic emission. The hatched area represents the pa- 
rameter space accessed by the diffuse flux constraints but 
not yet accessible by the point flux limits from direct 
searches. Therefore, diffuse flux limits can constrain the 
physics mechanism of neutrino production from individ- 
ual sources either to lower number density or to smaller 
fraction of power output of neutrino sources. The derived 
limits on the number density Ug of neutrino sources de- 
pend only on neutrino information, without making spe- 
ciflc associations with source class based on electromag- 
netic measurements. The region above the thick solid 
line is the excluded region by the upper bounds on the 
number density derived from the diffuse limits. 

The thick dark horizontal line in Fig. [T] and Fig. [2] 
indicates our primary constraint <i>c- We address the 
robustness of the constraint by focusing the discussion 
on the three assumptions involved in the calculation of 
Ng. 

The matter distribution within 5 Mpc of the Milky 
Way is far from uniform, which suggests the possibil- 
ity that the local number density of neutrino sources, 
n;, may be higher than the universal average of number 
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FIG. 3; The number density of neutrino sources Us plotted versus the expected neutrino luminosity predicted according to 
the fluxes of the model tested. The derived upper bounds from the diffuse flux show a stronger constraint than the limit from 
point flux from direct searches. The hatched area represents the limits accessed by the diffuse flux, but not yet accessible by 
direct measurement from the point source searches. Upper bounds on number density Us are computed for different neutrino 
sources (vertical arrows). Thin solid/dotted lines represent (diffuse (D) /point (P)) constraints on the number density with one 
and two orders of magnitude improvement. The region above the thick solid line is excluded by the diffuse flux limits. 



density, {Us). We argue that, in practice, the local inho- 
mogeneity affects only the class of sources characterised 
by low luminosities. The bright sources are too rare to be 
affected by local matter density variation - the likelihood 
of finding a bright neutrino source within 5 Mpc is small 
to begin with (if electromagnetic luminosity and neutrino 
luminosity are comparable), and the local enhancements 
in matter density insufhcient to change the probability 
of detection. 

On the other hand, if sources have a low mean lumi- 
nosity, then the nearest in the ensemble are more likely 
to be within a distance that could be affected by fluctu- 
ations in the local matter density. For example, within 4 
Mpc, the ratio between local matter density to the uni- 
versal average known as overdensity is estimated to be 
about 5.3 [30|. In this case, the flux constraint (Eq. [S]) 
should be adjusted to account for the higher density of 
local matter, $' = $ * (n//(ns))^/^. However, as Tab. HIl 
shows, the adjusted fluxes are below for a wide range 
of {L,y). For distances larger than 8 Mpc the overdensity 



of galaxies is rapidly approaching the universal mass den- 
sity. To exceed a source of a given luminosity L,^ must 
be within a distance di = (47r/3)^/^ • r^ax * {^'/^cY^'^. 
Assuming that the neutrino luminosity is comparable to 
the maximum luminosity in any electromagnetic band, 
no sources are found within a distance di that would vi- 
olate Ac- 



table IL Adjusted flux constraints "!>' to account for local 
enhancement of source density. 



erg/s 
6 X 10"^ 



GeV/cm^s 



d, 
Mpc 



GeV/cm-^s [30| 

4 X 10"'" SX 3 1.2 X 10"'-' 4 

2.5 X 10'*^ 7 X 10~'° 1.3 1.2 8.4 x 10'"' 8 

We address now the assumption that the neutrino lu- 
minosity distribution is proportional to a (possibly bro- 
ken) power law, which is observed for several classes of 
sources in the electromagnetic band. It was shown in 
[T^ l that Ng computed from the full distribution agrees 
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to within few percent with a simpler calculation using 
only the mean luminosity of the distribution. The rea- 
son is that the most common luminosities in the distribu- 
tion can only be observed at relatively short distances, so 
source evolution and cosmological effects are negligible. 
Sources with large luminosities are too rare to contribute 
significantly. On the other hand, it could be argued that 
the unknown luminosity distribution function is not well 
described by a (possibly broken) power law that typifies 
electromagnetic sources [IH. In this scenario, by using 
the limit on the maximum power density in Eq. 01 it 
is possible to constrain the mean luminosity for a given 
source class, if the number density is known, using the 
relation [T^ : 

< 47r — ln(— ) ^ = erg/s (9) 

For AGN selected in the x-ray band, Ug ~ 1.4 x 
10'' Gpc~^ [l^l, and the mean neutrino luminosity is 
< 2.4 X 10^^ erg/s, which is approximately two or- 
ders of magnitude lower than the average luminosity in 
the x-ray band. 

The constraint can be extended to energy spectra 
that differ from the assumed E^^ dependence, but the 
constraint applies over a restricted energy interval that 
matches the energy interval of the diffuse neutrino lim- 
its. Experimental diffuse limits span two different energy 
regions, VHE and UHE, and either limit can be inserted 
into Eq. [S] The restriction in energy range is required to 
avoid extrapolating the energy spectrum to unphysical 
values. In other words, for power law indices far from 
2, the spectrum must cut-off at high energies for indices 
7 < 2, or at low energies for indices 7 > 2. Subject 
to this restriction, we find that the constraint depends 
weakly on the assumed spectral index. For example, the 
constraints improve by a factor 2 for hard spectra (7 = 1) 
and weaken by roughly the same factor for soft spectra 
(7 = 3)[il. ^ 

To summarize, we have presented in this paper the 
constraint on neutrino fluxes from extragalactic point 
sources, which is E2(dN^/dE) < 1.4 x 10"^ (L^/2 x 
lO'*^ erg/s)i/3 GeV cm s . These constraints are 
one order of magnitude below current experimental lim- 
its from direct searches if the average distribution is 



comparable to the electromagnetic luminosity that char- 
acterizes the brightest AGN. As experimental data im- 
proves the derived constraints on fluxes from extragalac- 
tic sources E^(dNi^/dE) (x K'^(ff improves with the dif- 
fuse flux limits to the 2/3 power, while constraints on 
the number density oc Kdiff and the total neutrino 
power density V'^ oc Kdi / / improve linearly with the dif- 
fuse limits. We tested a number of model predictions for 
^'-point fluxes, and models which predict fluxes higher 
than the benchmark constraint have been restricted by 
this analysis. The constraint is strengthened for less lu- 
minous sources, and noncompetitive with direct searches 
for highly luminous explosive sources, such as GRBs. We 
found that the constraint is robust when accounting for 
the non-uniform distribution of matter that surrounds 
our galaxy, or considering energy spectra that deviate 
from -B^^, or various models of cosmological evolution. 
We also derived an upper limit on the maximum neutrino 
power density which is significantly below the observed 
power density from extragalactic cosmic rays. We showed 
that diffuse flux limits can strongly constrain the num- 
ber density of neutrino sources 11^. The constraints de- 
rived from the diffuse limits for sources with luminosities 
< 10^^ erg/s is stronger by few orders of magnitude 
compared to the point flux limits from direct searches. 
The parameter space accessed by the rig constrained from 
the diffuse limits for sources within this luminosity range 
is a challenge for direct point searches even for kilometer- 
cube neutrino detectors. The constraint suggests that 
the observation of extragalactic neutrino sources will be 
a challenge for kilometer scale detectors unless the source 
is much closer than the characteristic distance between 
sources, d/, after accounting for local enhancement of the 
matter density. Although the constraint cannot rule out 
the existence of a unique, nearby extragalactic neutrino 
sources, we note that assuming ~ L^, we found no 
counterparts in the electromagnetic band with the re- 
quired luminosity and distance to violate the constraint. 
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